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We show that the next-to-leading order perturbative prediction, matched with the next-to-leading logarithmic 
approximation for predicting both two-, three- and four-jet rates using the Durham jet-clustering algorithm, in 
the 0.001 < j/cut < 0.1 range gives a very accurate description of the data obtained at the Large Electron Positron 
Collider. This information can be utilized either for simultaneous measurement of the strong coupling and the 
QCD colour charges, or for improving the QCD background prediction in new particle searches at LEP2. 


1. INTRODUCTION 

Jet production rates provide one of the most in¬ 
tuitive ways to study the underlying parton struc¬ 
ture of hadronic events. However, lacking the 
necessary theoretical accuracy, for measuring the 
strong coupling a s , only the differential 2/3 distri¬ 
butions were used so far, where 2/3 is the event 
shape variable measuring the value of the jet res¬ 
olution parameter y cut for which the jet multi¬ 
plicity of a given event changes from three-jet to 
two-jet. The new next-to-leading order results for 
four-jet rates ]l] ||] raise the possibility of using 
LEP four-jet data for a s precision measurements. 

In electron-positron annihilation the widely 
known Durham j4|] jet-clustering algorithm has 
become an indispensable tool for classifying mul¬ 
tihadron final states into jets. This jet algorithm 
has the advantages of relatively small hadroniza- 
tion corrections and of the possibility of resum¬ 
ming large logarithms near the edge of the phase 
space (small y cut region), thus extending the va¬ 
lidity of the perturbative prediction. Recently a 
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new, Cambridge, jet-algorithm was proposed; it 
has similar resummation properties, but smaller 
hadronization corrections for mean jet multiplic¬ 
ities |fi. More detailed studies showed, however, 
that the small hadronization corrections found 
for the Cambridge algorithm in the study of the 
mean jet rate are due to cancellations among cor¬ 
rections for the individual jet production rates. 
Apart from the very small values of the resolu¬ 
tion parameter, y cut < 10~ 3 ' 2 , the Durham clus¬ 
tering shows, for the individual rates, comparably 
small (for y cut > 10 -2 ), or even much smaller, 
hadronization corrections |fi. We will be using 
data with y cut > 10 -3 ; therefore, in this talk we 
consider only multijet rates obtained using the 
Durham algorithm. 


2. THE THEORETICAL DESCRIPTION 


The n-jet rates are defined as the ratio of the 
n-jet cross section to the total hadronic cross sec¬ 
tion, and at next-to-leading order these take the 
general form 


Rn = 


O n— jet 
^tot 
n— 1 


+vM n ~ 1 0 


= »?(A) n 2 Bn(y C nt) (1) 

- 2)B„(y cut )Po ln(a; /i ) + C n (y cut )\ ■ 


In this equation ^(/i) = a s (/i)Cp’/27r, x M = 
y/y/s, where y is the renormalization scale and 
,/s is the total c.m. energy. The functions B n 
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and C n are independent of the renormalization 
scale, B n is the Born approximation and C n is 
the radiative correction. These functions in the 
case of the Durham clustering for n = 2 and 3 
were calculated based upon the ERT matrix el¬ 
ements [7] and for n = 4 they were obtained in 
Refs, fp] based upon the one-loop matrix ele¬ 
ments of Refs. H and tree-level matrix elements 
of Ref. (|). We use the two-loop expression for 
the running coupling, 


v(») = 


( 2 ) 
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{x = C A /C F and y t = T R N f /C F = N t /2C F ). All 
theoretical predictions in this contribution were 
obtained for five light-quark flavours at the Z° 
peak with M z = 91.187 GeV, T z = 2.49 GeV, 
sin 2 9w = 0.23 and a s (M z ) = 0.118. 

Multijet fractions decrease very rapidly with 
increasing resolution parameter y cut . Conse¬ 
quently, most of the available multijet data 
are at small y cut . It is well known that for 
small values of y cut the fixed order perturba¬ 
tive prediction is not reliable, because the expan¬ 
sion parameter (oc s /2tt) In 2 y cut logarithmically 
enhances the higher-order corrections. For in¬ 
stance, (a s /27r) In 2 0.01 ~ 0.4. Thus, one has to 
perform the all-order resummation of the lead¬ 
ing and next-to-leading logarithmic (NLL) con¬ 
tributions. This resummation is possible for the 
Durham algorithm using the coherent branching 
formalism Q. The two-, three- and four-jet rates 
in the NLL approximation are given in terms of 
the NLL emission probabilities Tj(Q, q) Q which 
have the following form: 

r q (Q,q)= (6) 
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We relate the rj(q) strong coupling appearing in 
the emission probabilities to the strong coupling 
at the relevant renormalization scale, ac¬ 

cording to the one-loop formula 

»?(P) 


v(q) = 


(9) 


where w(q, qo) was defined in Eq. (f|), and we use 
Eq. (j|) for expressing in terms of rj(M z ). We 
could also use a two-loop formula for 77 (g), but the 
result would differ only in subleading logarithms. 
However, we take into account a certain part of 
subleading soft logarithms with the inclusion of 
the K term. The K coefficient is renormalization- 
scheme dependent. In the MS scheme it is given 
by 0 


K = 
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The result of this resummation together with 
its renormalization-scale dependence in the case 
of four-jet rates was studied in Ref. where we 
found that the fixed-order and the NLL approxi¬ 
mations differ significantly. One expects that for 
large values of y cu t the former, and for small val¬ 
ues of ?/cut the latter is the reliable description; 
therefore, the two results have to be matched. 

The Durham multijet rates can be resummed 
at leading and NLL order, but they do not sat¬ 
isfy a simple exponentiation 0 (except for the 
two-jet rate). For an observable that does not 
exponentiate, the viable matching scheme is the 
R-matching ||, which we use according to the 
following formula: 


R 


R—match 


= R 


NLL 


( 11 ) 


y n - 2 (B n -B NLL 




c n - a 


NLL 


where B^ LL and C^ LL are the coefficients in the 
expansion of i?^ LL as in Eq. ([I]). 
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3. RESULTS 

In Figs. 1, 2 and 3, we show the theoretical 
prediction at the various levels of approximation: 
in fixed-order perturbation theory at the Born 
level (LO), at next-to-leading order (NLO), re¬ 
summed and R-matched prediction (NLO+NLL), 
and improved resummed and R-matched predic¬ 
tion (NLO+NLL+K), for the two-, three- and 
four-jet rates respectively. Also shown are the 
multijet rates measured by the ALEPH collabo¬ 
ration at the Z° peak 0 corrected to the parton 
level using the PYTHIA Monte Carlo [ p~2| ]. We 
used bin-by-bin correction and the consistency of 
the correction was checked by using the HERWIG 
Monte Carlo & The two programs gave the 
same correction factors within statistical errors. 
The errors on the data are the scaled statisti¬ 
cal errors of the published hadron level data, and 
we did not include any systematic experimental 
error, or the error due to the lradron-to-parton 
correction. In the lower parts of the plots we 
show the relative difference (data-theory) /theory, 
where theory means the ‘NLO+NLL+K’ predic¬ 
tion and we also indicated the renormalization 
scale dependence. 

Figures 1-3 deserve several remarks. First of 
all, we see that the inclusion of the radiative cor¬ 
rections improves the fixed-order description of 
the data, using the natural scale 2 + = 1 for larger 
values of y cu t- Secondly, the importance of resum¬ 
mation in the small y cu t region is clearly seen, but 
it is still not sufficient to describe the data at the 
natural scale, neglected subleading terms are still 
important. On the other hand, the improved re¬ 
summation seems to take into account just the 
right amount of subleading terms; this makes the 
agreement between data and theory almost per¬ 
fect over the whole y cu t region, as can be seen 
from the lower part of the plots. (In the case of 
four-jet rate, for y cut > 10 -1 ' 7 , 64 falls outside the 
±3 % band, one should keep in mind that in this 
region i) the renormalization scale dependence is 
relatively large and ii) the number of events is 
very small; therefore the statistical errors of the 
data and that of the hadron-to-parton correction 
are very large.) 

We found remarkably small scale dependence 



Figure 1. The QCD predictions for two- 
jet rates, with renormalization scale is set to 
2 + = 1, compared to ALEPH data. The lower 
part of the plot shows the relative difference 
S = (data-theory)/theory, where theory means the 
improved NLL approximation matched with the 
NLO result as explained in the text. The bands 
indicate the theoretical uncertainty due to the 
variation of the renormalization scale 2 + between 
0.5 and 2. 
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Figure 3. Same as Fig. 1, but for four-jet rates. 


for the ‘NLO+NLL+K’ predictions in the region 
y cut > ICC 3 . This feature, however, should be 
taken with care. On the one hand at any arti¬ 
ficial narrowing of the scale-dependence bands, 
e.g. at a crossover point, the bands almost cer¬ 
tainly do not represent the size of the trunca¬ 
tion error at that point. On the other, the im¬ 
provement, obtained by including the two-loop 
coefficient K, affects NNLL terms, but there are 
other contributions of the same order that are 
not taken into account (e.g. next-to-leading or¬ 
der running of a s and other dynamical effects). 
The scale dependence of the ‘NLO+NLL+K’ re¬ 
sult would consistently be under control only af¬ 
ter the inclusion of the complete set of NNLL 
terms. One expects however, that the scale de¬ 
pendence of the ‘NLO+NLL’ prediction is an up¬ 
per bound for the scale dependence of the per¬ 
turbative prediction, with subleading logarithms 
taken into account completely. Therefore, we 
may use our ‘NLO+NLL+K’ prediction for QCD 
tests, for instance, for measuring the strong cou¬ 
pling a s with the condition that we estimate 
the systematic theoretical uncertainty due to the 
scale dependence from the scale dependence of 
the ‘NLO+NLL’ prediction (obtained from vary¬ 


ing the scale between 0.5 and 2 as standard 
choice). The result of such a fit is given in Table 
1, where the central value was obtained using the 
‘NLO+NLL+K’ result with = 1 and the error 
represents 

• the statistical error on the data, 

• the systematic error due to changing the fit 
range (the whole range shown in Figs. 2 and 
3 was chosen) by one bin at both ends in 
both directions, 

• the error due to the use of different Monte 
Carlo programs (PYTHIA and HERWIG) 
for calculating the hadronization correc¬ 
tions, 

• and the error due to the variation of renor¬ 
malization scale as described above, 

all added in quadrature. This error is strongly 
dominated by the scale uncertainty. Of course, 
we could not include the systematic experimen¬ 
tal error. Also, in Table 1 we show the result 
of the fit when the renormalization scale is left 
as a free parameter. It is remarkable that, in 
the case of the three-jet rate, the natural scale 
Xfj = 1 is very close to the scale giving the small¬ 
est x 2 /d- 0 .f. (x^ = 0.92). On the other hand, in 
the case of the four-jet rate, the fitted scale is still 
somewhat lower than the natural scale. In our in¬ 
terpretation, this is due to the importance of the 
still neglected subleading terms. Indeed, we could 
also fit the four-jet data with the ‘NLO+NLL’ 
prediction, but with a very low scale ( 2 + ~ 0.2) 
and a slightly higher value of the strong coupling 
MAf z ) = 6.121). 

4. CONCLUSION 

In this talk we studied the perturbative descrip¬ 
tion of two-, three- and four-jet rates produced at 
LEP, obtained using the Durham clustering algo¬ 
rithm. We found that the best theoretical approx¬ 
imation that is currently available gives a remark¬ 
ably precise account of the data. In a previous 
publication Q, we found that the angular cor¬ 
relations defined on four-jet events are also well 
described (within ±5%) by perturbative QCD. 
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Table 1 

Results of a s (Mz) fits with fixed and fitted renor¬ 
malization scale using multijet data obtained by 
the ALEPH collaboration at LEP1. 



r 3 


& i?4 

a s (M z ) 

0.116 

0.1182 

0.1175 

X/j, fixed 

1 

1 

1 

X 2 /d.o.f. 

5.9/17 

29.2/12 

54.4/29 

errors 

statistical 

±0.0004 

±0.0003 

±0.0002 

fit range 

±0.0002 

±0.0002 

±0.0001 

hadroniz. 

±0.0008 

±0.0005 

±0.0001 

ren. scale 

±0.0015 

±0.0022 

±0.0017 

total 

±0.0018 

±0.0023 

±0.0018 

a s (M z ) 

0.116 

0.1176 

0.1173 

xn fitted 

0.92 

0.64 

0.7 

X 2 /d.o.f. 

5.5/17 

10.4/12 

29.6/29 


These observations suggest that the same proce¬ 
dure should provide an accurate prediction of the 
multijet backgrounds encountered in new particle 
searches and W-mass measurements at LEP2. 

We also performed a measurement of the strong 
coupling based upon multijet rates. We found 
a s (Mz) = 0.1173 ± 0.0018, where the error in¬ 
cludes the statistical and theoretical ones, but not 
the systematic experimental ones. 

These results were produced in part by a par- 
tonic Monte Carlo program called DEBRECEN 
, which is based upon the dipole formalism 
and can be used for the calculation of QCD 
radiative corrections to the differential cross sec¬ 
tion of any kind of infrared-safe three- and four- 
jet observable in electron-positron annihilation. 
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